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ABSTRACT 

Turbulence driven by magnetorotational instability (MRI) affects planetesimal formation by inducing dif- 
fusion and collisional fragmentation of dust particles. We examine conditions preferred for planetesimal for- 
mation in MRI-inactive "dead zones" using an analytic dead-zone model based on our recent resistive MHD 
simulations. We argue that successful planetesimal formation requires not only a sufficiently large dead zone 
(which can be produced by tiny dust grains) but also a sufficiently small net vertical magnetic flux (NVF). Al- 
though often ignored, the latter condition is indeed important since the NVF strength determines the saturation 
level of turbulence in MRI-active layers. We show that direct collisional formation of icy planetesimal across 
the fragmentation barrier is possible when the NVF strength is lower than 10 mG (for the minimum-mass solar 
nebula model). Formation of rocky planetesimals via the secular gravitational instability is also possible within 
a similar range of the NVF strength. Our results indicate that the fate of planet formation largely depends on 
how the NVF is radially transported in the initial disk formation and subsequent disk accretion processes. 
Subject headings: dust, extinction — magnetohydrodynamics (MHD) — planets and satellites: formation — 
protoplanetary disks — turbulence 



1. INTRODUCTION 



Formation of kilometer-sized planetesimals is the initial 
step of planet formation in protoplanetary disks. Several 
mechanisms have been proposed for planetesimal formation, 
whic h include gravitational in s tability (GI) of dust s ubdisks 
(e.g.. iGoldreich & Ward| [l97l Uohansen et alj|2007t lYoudinl 



1993 



201 1 ) and direct coa gulation (e.g., Weidenschilli ng & Cuzzil 
lOkuzumi et alJl20T2h . However, the outcome of these 



processes strongly depends on the turbulent state of the gas 
disk. Turbulence is know to concentrate particles of par- 
ticular sizes, which could assist their g ravitational collapse 
dCuzzi et alJl200U Uohansen et al.l l2007). On the other hand, 
turbulence also stirs up dust subdisks and thereby stabilizes GI 
dTurner et alJl20ldh . In addition, the relative velocity induced 
by turbul ence can lead to catas trophic disruption of large dust 
particles ( Uohansen et alj|2008b . 

The magnetorotational instability (MRI; Balbu s~& Hawlevl 
1991) has been thought as the most plausible driv- 
ing mechanism of protoplanetary disk turbulence. The 
activity of MRI largely dep ends on non-ideal MHD 
effects (e.g.. ISano & Mivamal [19991 iBai & Stond IMlt 
iWardle & Salmeronll2012l) A high ohmic resistivity creates 
an M RI-inactive "dead zone" near the midplane (Gammie 
1996), which reduces turbulence strength. Importantly, the 
size of the dead zone depends on the amount of tiny dust parti- 
cles, because they effectively red uce the gas ioniza tion degree 
and hence enhance the resistivity (ISano et al.l2000l) . Although 
the dead zone has often bee n invoked as a favorable site for 
planetesimal formation (e.g.. lCieslall2007HBai & Stondl2010t 
lYoudinll201 ll) . self-consistent modeling of the MRI-dust co- 
evolution has not been done so far. 

Another important, but much less appreciated factor is the 
net vertical flux (NVF) of the magnetic fields. The NVF 
indeed matters since it determines the saturation level of 



MRI- driven turbulence (e.g jHawlev et alj!9 95: Suzu ki etaT] 
2010). This is especially true when a large dead zone is 
present, for which case the vertically integrated accretion 
stress decreases with de creasing N VF dGressel et all 120121 
lOkuzumi & Hirosdl201l[ henceforth lOH 1 11) . 

In this Letter, we investigate possible pathways of planetes- 
imal formation taking into account the dependence of MRI- 
driven turbulence strength on dust size distribution and NVF 
strength. In our previous paper (OH11), we systematically 
studied the saturated state of MRI-driven turbulence with lo- 
cal stratified, ohmic -resistive MHD simulations. We obtained 
an analytic prescription for the saturation level as a func- 
tion of the ohmic diffusivity and NVF strength. Using this 
prescription together with ionization balance calculation in- 
cluding grain charging, we are able to determine turbulence 
strength consistently with the amount of tiny grains and NVF 
strength. We consider the formation of rocky and icy plan- 
etesimals independently, because icy particle s have a high 
sticking efficiency compared to rocky particles (Choks hi et al.l 
119931: iGundlach etai1l201 ll) . We test direct collisional for- 
mation of icy planetesimals outside the snow line, and rocky 
planetesimal formation via GI inside the snow line. 

2. MODEL DESCRIPTION 

We consider a protoplanetary disk around a solar mass star. 
We take the gas surface densit y £ fl from the minimum-mass 
solar nebula (MMSN) model of Ha vashil (119811) . The gas tem- 
pe rature is taken from the pass ive, optically thick disk model 
of IChiang & Goldreich (1997). This model well approxi- 
mates the gas temperature inside a dead zone, since MRI- 
driven heating occurs mainly at the upper bo undary of the ac- 
tive l ayer where the optical depth is small dHirose & Turner 
1201 ll) . The assumed temperature gives a snow line at orbital 
radius r ~ 1 AU. The gas density p q depends on the distance 
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z from the midplane as p g - (L g / ^Jlnh) exp(-z 2 /2/i 2 ), where 
h is the gas scale height given as the sound speed c s divided 
by the Keplerian frequency D.. Because of gas pressure sup- 
port, the gas disk rotates at a slightly sub-Keplerian velocity. 
This causes systematic radial drift of dust particles relative to 
the mean gas motion. The drift speed reaches the maximum 
|frfr,maxl ~ 30 m s _1 when the stopping time t s of the parti- 
cle equals Q _1 dWeidenschillindl 19771) . At r ~ 1-5 AU, the 
dimensionless stopping time Of s = 1 corresponds to particle 
radius a ~ 1 m. 

We consider MRI-driven with a dead zone. The most im- 
portant turbulent quantity for dust evolution is the random ve- 
locity of the gas since it determines the collision velocity and 
diffusion coefficient of dust particles. IOH1 ll found that the 
gas velocity dispersion at the midplane, <5% m id, is well ap- 
proximated as 



a = 0.1 fim 



Sv, 



'#,mid 



= -\/0/78acoreC s 



(1) 



where a core is the accretion stress integrated over low altitudes 
and normalized by £ 9 c 2 . Equation holds no matter if the 
midplane is magnetically dead, since hydrodynamical waves 
created in active layers propagate across dead zone bound- 
aries. lOHl ll also found that ff cole is determined by the dead 
zone size and NVF as 



510 / 0.54/i res \ / 3.6/ia 
exp +0.011 exp 



(2) 



where ySNVF.mid = 87rp ?jlll idC 2 /B 2 vf * s tne midplane plasma 
beta defined by the NVF strength Bnvf, and h KS and are 
quantities that characterize the vertical extent of the dead 
zone. Precisely, h KS is the height below which the char- 
acteristic MRI wavelength A KS = 2nrj/vA Z in the resistive 
MHD limit exceeds h, whereas h\ is the height below which 
the ohmic Elsasser number A = v\ JrjQ. falls below unity, 

where va z = Bnvf/ y/^npg is the Alfven speed defined by the 
NV F and 77 is the ohmic diffusivity. Linear stability analy- 
sis dSano & Miv ama 1999) shows that ohmic resistivity sup- 
presses the most unstable MRI mode at A < 1 (z < h^) but 
less unstable modes survive as long as A res < h (z > h KS ). 
Thus, the region /i res < z < h\ can be interpreted as the tran- 
sition layer between the active and dead zones. The active 
layer extends up to z = &ideai> above which MRI is stabilized 
becaus e of a low local plasma beta dSano & Mivamal [19991 
lOHl ll) . In the saturated state, /zideai is approximately given by 
fy dfiai = \2 ln(yS N vF,mid/3000)] l/2 h according to the simulations 



Iidj 

of lOHTll 

For given (5u ff , m id, the gas random velocity 5v g at each height 
z is given as 6v g = exp(z 2 /4/i 2 )c)u 9jln id (lOHl lh . The vertical 
diffusion coefficient for du st is given by D z = D z o/[l +(Qf 5 ) 2 ] 
dYoudin & Lithwickll2007l) . where 



D z{) * 6v z g /3Q 



(3) 



is the verti cal diffusion coefficient for passive (Q. t s <K 1) con- 
taminants (Fromang & Papaloizou 200(J lOHl lh . Note that 
Equation (01 applies even in dead zones, which is consistent 
with the fact that hydrodynamical waves propagati ng from ac- 
tive la yers have a finite correlation time ~ £2 _1 dGressel et al.l 
120121) . The turbulence-driven relative velocit y Av t of dust par- 
ticles is calculated from the prescription of lOrmel & Cuzzil 
d2007l) for Kolmo gorov turbulence. In reality, a Kolmogorov 
energy cascade may not be established for waves in dead 
zones. However, the assumed Au, at least gives a reasonable 
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Figure 1. Example of the dead zone size. The dot-dashed curve shows z = 
''ideal for Bnvf = 5 mG, while the solid and dashed curves show z = h\ and 
Z = h ms , respectively, for Bnvf = 5 mG, Z sma u = 10~ 4 , and a = O.lfrni. For 
reference, z = h and 3/i are shown by the dotted curves. 

estimate of the relative velocity for marginally or fully decou- 
pled (Qf s > 1) particles since their relative velocity is deter- 
mined by t he largest-scale gas m otion with correlation time 
~ Q" 1 (see lOrmel & Cuzzill2007l) . 

Equations (Q]) and (f2j) predict <Su 9 , m id as a function of h les , 
/ia, and Bnvf- We calculate the critical heights by con- 
sidering the ionization balance at each z taking into ac- 
count grain charging. We use the analyt ic solution of th e 
ionization balance equations derived by lOkuzumil (I2009I) . 
which gives the ionization degree for arbitrary dust size 
distribution. T he ionizing sources we inclu de are Galac- 
tic cosmic rays (lUm ebavashi & Nakano 2009), stellar X-rays 
(llgea & Glassgoldlll999l: iBai & Goodman! I2009I) . stellar en- 
ergetic protons (|Turner & Drakel I2009I) . and radionuclides 
(Umebayashi & Nakano 2009). Inclusion of cosmic and stel- 
lar protons gives the minimum estimate of the dead zone size 
since strong T-Tauri winds may shield these particles well 
above disk surfaces. From th e ionization balance, we calcu- 
late the vertical profile of 77 (Blaes & Balbus 1994), and ob- 
tain hies an d f° r given Bnvf- Figure Q] plots z - /i re s an d 
Z = h\ versus r for Bnvf = 5 mG assuming that 0. 1 fim sized 
grains are uniformly mixed in the gas with mass abundance 
Zsmaii = £ sma u/E 9 = 10~ 4 . Note that the dead zone size de- 
pends on Bnvf; the larger Bnvf is, the smaller /i res and /ia are. 

3. ICY PLANETESIMAL FORMATION ACROSS THE 
FRAGMENTATION BARRIER 

Planetesimal formation via direct coagulation is limited by 
the fact that marginally decoupled (Qf s ~ 1, a ~ 1 m) dust 
aggregates experience high-speed collisions that can lead to 
catastrophic disr uption (the so-called fragmentation barrier; 
iBrauer et al.H2.Q08D . If the disk is laminar, the maximum col- 
lision velocity is determined b y the differential radial drift 
velocity ~ |^ r , ma xl ~ 30ms~' (Braueretal. 2008). Recent 
numerical collision experiments (IWada et al.ll2009l) show that 
aggregates made of submicron-sized icy grains stick at col- 
lision velocities up to « 50 m s -1 , suggesting that icy plan- 
etesimal formation via direct coagulation is possible in lami- 
nar disks. However, if the disk has MRI-active layers, a ran- 
dom velocity of ~ £i> o,mid is added to the co llision velocity for 
Qf s ~ 1 aggregates (Ormel & Cuzzi 2007). The question is: 
can icy aggregates grow across the fragmentation barrier even 
if the MRI-driven turbulence enhances the collision velocity? 

To get a feeling of how the fragmentation barrier depends 
on dust size distribution and NVF, we begin with a simple 
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Figure 2. Midplane gas random velocity (JlVmid at 5 AU vs. the mass abun- 
dance Z smii n of small grains for different values of the NVF strength Bnvf- 
The filled circles and stars mark Z sraa n below which h reli and h\ vanish, re- 
spectively. Above the thick gray line, the total collision velocity for large 
(Hf , l an ,e ~ 1) aggre gates exceeds the catastrophic disruption threshold for 
ice {Wada et al. 20Q3). 

two-population model in which large, marginally decoupled 
(Qf. S; i alg e ~ 1) aggregates coexist with 0.1 fim sized small 
grains. We calculate 6v g ^ assuming that only the small 
grains contribute to the ionization balance. The mass abun- 
dance Z srrm ii of the small grains is taken as a free parameter. 
Figure |2] shows (5i> ff , m id as a function of Z sma u for different val- 
ues of Z?nvf- For fixed Bnvf, (^.mid increases with decreasing 
Zsmaii because the dead zone is smaller when small grains are 
less abundant. The thick gray line shows tU> ff , m id = 38 m s , 
above this line, the total collision velocity for two Q.t s ~ 1 
aggregates, Aui al - ge ~ (Su 2 g mid + f^. max ) 1/2 , exceeds the disrup- 
tion v elocity i>disr = 50 m s _1 for icy aggregates (Wa da et al.1 
2009). If Bnvf — 1-1 mG, Aui arge falls below Udisr for all val- 
ues of Z sma n. The reason is two-fold. For such small Z?nvf, a 
substantially large (h^ ~ h) dead zone is present even in the 
absence of small grains (note that /ia is larger when Z?nvf is 
smaller). In addition, the small Bnvf leads to a low saturation 
level in the upper active layer, resulting in a low gas veloc- 
ity dispersion (note that <5t^. m id °c B NV f when /ia » h/4). 
However, if Bnvf goes up to 3.2 mG, Ai> lal - ge exceeds Udisr for 
Z sma u < 10~ 5 , for which the active layer reaches the midplane 
(/za = 0). For Bnvf S 10 mG, Aui arge exceeds Ddisr even if 
Z sma ii is as large as the interstellar value 10~ 2 . 

For Bnvf ~ 3 mG, large aggregates can grow across 
the fragmentation barrier only if small grains are sufficiently 
abundant. To show that it is indeed possible, we simu- 
late the evolution of full dust size distribution at 5 AU us- 
ing a coagulation-advection-diffusion equation (e.g., ICieslal 
2007). In this simulation, we begin with 0.1 fim icy dust 
grains well mixed in the vertical direction, and follow the 
evolution of the vertical dust size distribution due to co- 
agulation/fragmentation and vertical settling/diffusion. The 
collision velocity Ay of dust particles takes into account 
Brownian motion, sy stematic drift in the mean gas motion 
dNakagawa et al.l 119861). an d turbulence-driven random mo- 
tion (lOrmel & Cuzzil 120071) that depends on Sv g . We take 
6v g to be time-dependent, namely, to be consistent with the 



dead zone size c alculated from the full dust size distribution 
(Okuzumi 2009) at each time step. The total dust mass abun- 
dance Z = is taken to be 10~ 2 , and the loss of dust 
materials due to the radial drift is neglected. Indeed, coagu- 
lation proceeds faster than radi al migration if the dus t grows 
into highly porous aggregates (Okuzu mi et al.l 120121) . How- 
ever, in order to make our model as simple as possible, we 
neglect both porosity evolution and radial migration. The 
mass of an aggregate after collision is given by m, + sm p , 
where m, and m p (< m,) are the masses of the target and pro- 
jectile, and s(< 1) is the dimensionless sticking efficiency 
that depends on Av. We assume 5=1 for Au < 10 m s 
and i=l- ln(Au/10 ms _1 )/ln5 for Ay > 10 ms~' in ac- 
cordance with the r esult of numerica l collision experiments 
for icy aggregates (Wad a et all 120091) . The target loses its 
mass (s < 0) when Ay is larger than the disruption threshold 
fdisr = 50ms _1 . The fragments (whose total mass is (l-s)m p ) 
are assumed to be in the form of 0.1 fim sized constituent 
grains. This assumption is quite simplistic, but still is reason- 
able as a first-order approximation since the total mass of frag- 
ments in aggregate-a ggregate collision tends to be domina ted 
by the smallest ones (Wurm et al. 2005; Wada et al. 2009). 

Figure [3] shows the simulation result for Bnvf = 3.2 mG. 
The top panels display the vertical dust mass density pd per 
unit logarithmic particle radius log a at different times t after 
the initial state. Initially, dust particles grow without catas- 
trophic disruption because of low collision velocities. As the 
particles grow, disruption becomes more and more significant, 
because the collision velocity increases with increasing stop- 
ping time and because the dead zone shrinks with the deple- 
tion of small grains. The second effect becomes prominent 
when t x 560 yr, at which the dead zone disappears and 
consequently 6v g ^ m i,i reaches 55 ms _1 . However, the large 
amount of tiny grains produced by the catastrophic disrup- 
tion are quickly diffused to high altitudes, reduce the ioniza- 
tion degree, and finally "revive" the dead zone. The revived 
dead zone suppresses £>y ff>m id down to 27 m s~', which is low 
enough for larger aggregates to continue growing. In this way, 
larger aggregates grow beyond the meter-size fragmentation 
barrier. 

The middle panels of Figure|3]plot the dust mass abundance 
Z per unit log a, while the bottom panels show the vertical ge- 
ometric optical depth T geo of dust per unit log a. We see that 
the dust mass is dominated by large aggregates while the op- 
tical depth is dominated by small fragments. This indicates 
that small fragments carry a minor fraction of the total dust 
mass but still help the grow of larger aggregates by providing 
a large surface area needed to maintain a low enough ioniza- 
tion degree. 

4. ROCKY PLANETESIMAL FORMATION VIA SECULAR GI 

Coagulation of rocky aggregates is severely restricted by 
the fragmentation ba rrier since the disruption threshold i s 
as low as 1-5 ms _1 (IWada et al.ll2009t iGutfler et al.l l2010h . 
Using that the radial d rift speed \vd r \ is appr oximated as 
2[i>rf r ,maxl^fj for £lt s <k 1 ( Weidenschillingl l 1 9771) and assum- 
ing that Udisr ~ 5 m s~' and |tw,maxl ~ 30 ms _1 , we find 
that \vdr\ reaches i>disr when Qf, ss 0.08, which corresponds 
to a ~ 10 cm at ~ 1 AU. Hence, a simple coagulation sce- 
nario does not account for rocky planetesimal formation even 
without turbulence. 

One mechanism that can lead to rocky plane tesimal forma- 
tion is the so-called secular GI (Youdin 2011). It is a grav- 
itational collapse of dust materials driven by the combina- 
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Figure 3. Result of a coagulation simulation for icy dust at 5 AU for Bnvf = 3.2 mG. The top density plots show the dust mass density p r ; per unit logarithmic 
particle radius log a at different heights z and different times t. The solid and dashed lines in the top panels mark the dead zone critical heights z = h\ and z = h ms , 
respectively. The middle and bottom panels show the mass abundance Z and geometric optical depth r geo of dust per unit log a, respectively. Shown in the top of 
the middle panels is the dimensionless stopping time n?, m y of the dust particles at the midplane. 

tion of self-gravity and gas friction. An important feature 
of the secular GI is that it works even if Q.t s <K 1. Thus, 
the secular GI can allow gravitational collapse of dust par- 
ticles whose growth is limited by the fragmentation barrier. 
Instead, the secular GI requires sufficiently weak radial dust 
diffusion in order for the particles to collapse faster than they 
drift inward. For Qf s « 0.1, the radial diffusi on coefficient 
D r must be lower than 10~ 5 c s h at the midplane (Youdi nl201 It 
Take uchi & Ida1l2012h . 

To assess whether the secular GI operates for the 
fragmentation-limited aggregates, we estimate the MRI- 
driven diffusion coefficient assuming that the aggregates coex- 
ist with 0.1 j-im sized fragments of mass abundance Z sma n. We 
also assume that D r * D z , which is » D z q for Qf s » 0.1 <K 1 
(see Equation (O). Figure [4] shows the midplane radial dif- 
fusion coefficient Z) r m jd at 1 AU as a function of Z sma u for 
different values of Bnvf- We find that D ri mid exceeds 10 -5 c s h 
for Z?nvf ^10 mG even if Z sma u is as large as the interstellar 
value 10~ 2 . Thus, the secular GI of the fragmentation-limited 
aggregates requires Z?nvf ^ 10 mG. This requirement is sim- 
ilar to to that for direct icy planetesimal formation shown in 
Section[3] 
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5. CONCLUSION AND DISCUSSION 

We have investigated how planetesimal formation depends 
on the amount of tiny grains and the strength of the NVF For 
MMSN disks, we have shown that the existence of a large 
dead zone and an NVF weaker than 10 mG is preferable for 
planetesimal formation via both direct coagulation and secu- 
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Figure 4. Midplane radial diffusion coefficient Z) r m id at 1 AU vs. the mass 
abundance Z sraa n of small grains for different values of the NVF strength 
-Snvf- Below the thick gray line, the fragmentation-limited rocky aggregates 
(Clt s larg e ~ 0. 1) experience secular gravi tational collapse faster than they drift 
inward I Youdin 201 1: Takeuc hi & Idal2012h . 

lar GI. The obtained criterion for the NVF depends on the disk 
surface density S 9 , since 5v g m {& oc Bnvf/ ^f^g in the presence 
of a large dead zone (see Equations ([T} and (O). If E 9 is 10 
times larger than the MMSN value, then the upper limit on 
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Snvf goes up to 30 mG. 

We have neglected the effects of non-ohmic magnetic dif- 
fusivities. Ambipolar diffusion may stabilize MRI near the 
upper boundary of the active layer (Bai 201 1). The effect of 
Hall diffusion is more uncertain; it can stabilize or destabi- 
lize MRI dep ending on the sign of NVF relative to the disk 
rotation axis (IWardle & Salmeronll2012l) . Inclusion of these 
effects may change our results quantitatively, but the general 
trend that a weak NVF is preferable may be unchanged. 

Weak turbulence is also beneficial to the growth of solid 
objects larger than planetesimals. Density fluctuations in 
turbulence gravitationally interact with planetesimal-size ob- 
jects and thereby enhance their collision velocity. If MRI is 
fully active, the resulting gravitational stirrin g likely causes 
catastrophic disruption o f planetesi mals (e.g.. Ilda et al.l 20081 : 
iNelson & GressefcOlOh . However. iGressel et all (120 1 21) have 
recently shown that even weakly bound planetesimals are able 
to grow if a dead zone is present and if the NVF is weaker than 
3 mG (assuming the MMSN surface density). Thus, a weak 
NVF and a wide dead zone are preferable for the growth of 
solid bodies up to protoplanets. 

A realistic range of the NVF strength is poorly constrained 
by direct observations since the differential rotation of the disk 
can produce toroidal mag netic fields as strong as 0. 1-1 G even 
inside the dead zone dTurner& Sanol 12008). An important 
property of NVF is that the total NVF is a conserved quan- 
tity of a magnetized disk. Thus, the total NVF of a proto- 
planetary disk is directly determine d by how the disk formed 
from weakly magnetized (~ 10 /uG: iHeiles & Crutcherll2005l) 
molecular clouds. Nevertheless, what fraction of the magnetic 
flux is brought to the planet-forming inner disk region is not 
evident because non-idea l MHD processes a lso work during 
the disk formation (e.g., Machida et al.H201~Tb . Furthermore, 
on longer timescales, the NVF may be radially transported 
due to inward mass accretion (Rothstein & Lovelace 2008) 
and/o r outward macroscopic magnetic diffusion (Lubow et al. 
1994) in the active layers. So far, the origin and global trans- 
port of NVF has not been paid attention in the context of 
planet(esimal) formation. We hope that this Letter will en- 
courage further investigation on this issue. 
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